Three-body Coulomb breakup of ^^Li in the complex scaling method 
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Coulomb breakup strengths of '^^Li into a three-body ^Li+n+n system are studied in the complex 
scaling method. We decompose the transition strengths into the contributions from three-body 
resonances, two-body "^"Li-l-n" and three-body "®Li-|-n-|-n" continuum states. In the calculated 
results, we cannot find the dipole resonances with a sharp decay width in "Li. There is a low 
energy enhancement in the breakup strength, which is produced by both the two- and three-body 
continuum states. The enhancement given by the three-body continuum states is found to have 
a strong connection to the halo structure of ^^Li. The calculated breakup strength distribution is 
compared with the experimental data from MSU, RIKEN and GSI. 
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Studies of unstable nuclei have obtained much atten- 
tion with the development of radioactive beams T|. The 
^^Li nucleus is known as a typical Borromean system, in 
which the ^Li-fn and n+n subsystems do not have any 
bound states, but the total ^Li-|-n-|-n system has a bound 
state. This Borromean mechanism is considered to play 
an important role in the formation of a halo, but it is not 
yet fully understood. The observed large matter radius 
of "'^"'^Li, which is an evidence of the halo, suggests a large 
mixing of the (lsi/2)^ neutron component in addition to 
the one. Another interesting problem related to 

the halo structure of ^^Li is a characteristic property of 
the excitation mode. For the excited states of ^^Li, the 
so-called soft dipole resonance[2 is expected in the 
low-energy region. In the shell model picture, the ma- 
jor component of the soft dipole resonance is described 
as (lsi/2)(0pi/2)- Thus, the behavior of the Is- and Op- 
orbits of valence neutrons is very crucial to understand 
the halo structure and the excited states in ^^Li. 

Experimentally, measurements of the Coulomb 
breakup strength distributions of ^^Li have been carried 
out by three groups at MSU 4], RIKEN Q and GSl|l. 
The low energy enhancement of the strength seems to 
indicate the existence of the soft dipole resonance, al- 
though the shapes of distributions obtained by the three 
experiments at different incident energies do not coincide 
with each other. In addition to these measurements, 
observations of the two-body correlation provide us with 
a key to discuss the mechanism of the breakup reaction. 
The measured invariant mass spectrum of ^Li-|-n shows 
the low energy enhancement 0, Q- This result implies 
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the existence of the low-lying Is-orbit near the ^Li-|-n 
threshold energy in ^"Li. Such a Is-orbit in the ^Li-n 
system is expected to provide "a continuum structure" 
in the three-body breakup reaction of ^^Li. On the other 
hand, in the case of the two-body breakup reaction 
of ^^Be, the low energy enhancement observed in the 
strength is understood as "a continuum response" of a 
large low-momentum component in the halo structure 
of the ground state Q. In order to understand the 
observed enhancement of the ^^Li breakup strength, 
therefore, we must consider the continuum structure 
of the ^^lA-^-n binary component and the continuum 
response of ^Li+n+rt in addition to the three-body 
resonance. 

Theoretically, many methods such as the Faddeev 
method, the hyperspherical harmonics approach and so- 
phisticated variational methods have been developed to 
solve the Borromean systems 9]. However, there is a dis- 
crepancy between some theoretical results for the soft 
dipole resonances in ^^Li. Garrido et al.lO] calculated 
the dipole strength distribution, and predicted at least 
three dipole resonances, but we did not obtain any reso- 
nant solution with a sharp width (r/2 <Er) in the vari- 
ational methodfOi This discrepancy is considered 
to come from a difference on the ^Li-n potential. The 
dipole strength distribution calculated by Garrido et al. 
does not agree with the observed strengths within the ex- 
perimental ambiguity. The peak energy is slightly lower 
and the width is much sharper than the experimental 
one. Thus, it is strongly desirable to obtain detailed in- 
formation about the ®Li-n potential, the excited resonant 
states and the continuum responses through analyses of 
the Coulomb breakup reaction of ^^Li. 

For this purpose, we have been developing the appli- 
cability of the complex scaling method (CSM)jl3| which 
has recently received much attention for finding three- 
body resonances, 10,. 12j . It is a big advantage of GSM 
that for an unbound system conduces to the separation 
not only between resonances and continuum states but 
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also between different kinds of continuum states starting 
from different thresholds^]. This advantage of CSM is 
exploited in the calculation of transition strengths of un- 
bound states beyond the two-body systems. We showed 
a successful result for the three-body Coulomb bre akup 
reaction of a simpler two- neutron halo nucleus ^He[l4j. 
The results for ^He are summarized as (i) no resonance 
peak corresponding to the soft dipole mode is obtained, 
and (ii) the ^He(3/2~)-|-n binary component dominates 
the El strength distribution and the responses of the 
other components, such as three-body continuum states 
of ^He-l-n-l-n, are very small. 

In this letter, we extend this method to ^^Li and inves- 
tigate the continuum structures and responses through 
the Coulomb breakup reaction. We briefly explain an ex- 
tended ^Li+n+n model for ^^Li, and report the results of 
the Coulomb breakup strength distributions comparing 
with experimental data. From the obtained results, we 
discuss the mechanism of the breakup processes. 

We describe ^^Li with an extended ^Li-l-n+n three- 
body model 15]. The Hamiltonian of this model is given 
in the orthogonality condition model^l6] as follows: 

3 2 
if ("Li) = HCU)+Y,t^~TG + Y,Vcn.+Vnn 

1=1 i=l 

+ Ap_F \4'pf)(4'pf\, (1) 

where i?(^Li), ti and Tq are the internal Hamiltonian of 
^Li, the kinetic energy of each cluster and the center-of- 
mass of the three-body system, respectively. The ^Li-n 
potential, Vc„, is given by a folding-type one with MHN 
interactional?.]. For the potential for two valence 
neutrons, the Minnesota potential[ig is used with pa- 
rameter M=0.95. The last term in Eq. ([^I is a projection 
operator to remove the Pauli forbidden (PF) states from 
the ^Li-n relative motionll9|. The definition of the PF 
states is also given in Ref. (l5|, where the value of Xpp 
is taken as 10^ MeV in the present calculation. 

The folding-type ^Li-n potential was originally con- 
structed so as to produce energy splittings in the ^°Li 
spectra, such as 1+-2+ (for the pi/2-neutron) and l~-2~ 
(for the Si/2-neutron) due to the coupling between spins 
of the valence neutron and ^Li(3/2)^]. However, in the 
study of ^^Li0, we discussed to add a phenomenologi- 
cal tail potential to the original folding-type potential to 
improve the behaviour of the tail part of the ^hi-n poten- 
tial. The behavior of the s-wave state near the threshold 
is very sensitive to the tail part of the potential due to the 
spatial extension of the wave function. The tail poten- 
tial also plays an important role in lowering the energy 
of the (lsi/2)^ component with respect to that of the 
(Opi/2) ^-component in -'^-'^Li(^Li-l-n-l-n). Then, we have 
two parameters in the ^Li-n potential; one is the S to 
change the strength of the second range of the folding 
part, and another the strength of phenomenological tail 
potential given by a Yukawa- form][l5j. They are deter- 
mined so as to reproduce the l"*" resonance of ^°Li at 0.42 



MeV|20j and the s-wave property which is a virtual state 
showing a large negative scattering length. 
The wave function of ^^Li is given as 

*-'("Li) = 5]^{[cI>3/2-(C.),xl(nn)]^}. (2) 

i 

Here, the ^Li nucleus is expressed by a multi- 
configuration tti <i>'^/^ (Ci) in order to take into ac- 
count the neutron pairing correlation [22. Since the mix- 
ing amplitude depend on the relative distances be- 
tween ^Li and two valence neutrons in the ^Li-|-n-|-n sys- 
tem, we express it by the function x:^(nn). We should 
notice that xi (nn) — > x (plane wave) at large distances 
between ^Li and two valence neutrons. For the total spin 
J of ^^Li, j expresses the spin of two valence neutrons, 
and then J = 3/2 (g) j. 

We describe the wave function of valence neutrons us- 
ing the combined set of the two kinds of the basis states; 
the cluster orbital shell model (COSM; V-type) and the 
extended cluster model (ECM; T-type), which we call as 
the hybrid-TV modeljl (IB Hj: 

xKnn) = xly(lv) + xlT(^T), (3) 

where and are V-type and T-type coordinate sets, 
respectively. The radial part of the relative wave func- 
tion is expanded with a finite number of Gaussian basis 
functions centered at the origin. 

By employing Ci = (0^3/2)^^ and C2 = 
(0^3/2)^(0^1/2)^ for p-shell neutrons in ^Li {Nc is 
2), we solve a coupled-channel ^Li+n+n three-body 
problem. We use the MHN interaction^^] to calculate 
the neutron pairing correlation in ^Li, which leads to 
the ]a2]^ = 15 % of the pairing excitation J^. When 
the valence neutrons approach to ^Li, the coupling 
between the motion of the valence neutrons and the 
pairing correlation in ^Li becomes stronger. It was 
discussed that this dynamical coupling in ^*'Li(^Li+rt) 
provides the so-called pairing-blocking effect 22j which 
explains the lowering of Is-orbits in ^'^Li. For ^^Li, we 
adjust the (Op)^-(ls)^ pairing coupling between valence 
neutrons to reproduce the observed binding energy of 
"Li (O.31MeV)0. 

In the analysis of excited states in ^^Li, we prepare 
the three types of the ^^Li wave function; P-1, P-2 and 
P-3, which are characterized by the (lsi/2)^ probability 
in the ground state. In Table HI we list the properties 
of the prepared ^"^Li and ^"'^Li wave functions. The s- 
orbit properties are different among them, which affect 
on the scattering length of the "'^^Li s-state and the size 
of the halo structure of the "'^^Li ground state. We also 
pay attention to the effects of the halo structure on the 
breakup strength. They can be seen from the responses 
of resonances and continuum states in ^^Li. 

We calculate the unbound states of ^^Li applying CSM 
to the extended ^Li-f n-t-n system, where the relative co- 
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TABLE I: Results for the three types of wave functions of 
the present model; (upper) Scattering lengths and ener- 
gies E of the virtual state (2~,1~) of ^"Li, (lower) (lsi/2)^ 



probability P and matter radius R 


of the 


^^Li ground state. 


lOLi 


P-1 


P-2 


P-3 


a. (2-) [fm] 


-12.7 


-17.0 


-21.7 


E{2-) [MeV] 


-0.05 


-0.03 


-0.02 


a.(l-) [fm] 


-6.5 


-8.6 


-10.7 


E{1-) [MeV] 


-0.08 


-0.06 


-0.05 



'Li 



P[(lSl/2)'] [%] 

-Rm [fm] 



P-1 


P-2 


P-3 


21.0 


29.4 


38.8 


3.33 


3.58 


3.85 



exp. 



3.12±0.16" 
3.53±0.06'' 



" Reference |2J|, " Reference [ga] 

ordinates between ^Li and the two neutrons are trans- 
formed with a scaling angle 6 as 



,V.T 



lV.T< 



(4) 



The momenta corresponding to the asymptotic channel 
a of ^Li+n+n and ^*'Li'^*-'-|-7i are also transformed as 



(5) 



Here, we notice that 29 corresponds to a rotation an- 
gle of the cuts in the Riemann sheets of complex ener- 
gies, and the angular part of the wave function does not 
change in CSM. In Fig. ^ we show a schematic energy 
eigenvalue distribution of the complex-scaled ^Li+n+n 
system governed by the ABC-theorem^^. When 6 = 0, 
(unbound) scattering states are obtained on the real en- 
ergy axis, which includes all components of resonances 
and continuum states. For a finite value of 9, the contin- 
uum states are obtained on the Riemann cuts rotated 
down by 29, and hereafter we call these rotated con- 
tinuum states as the continuum ones. When we take 
a large 6, as shown in Fig. Q in addition to the three- 
body bound states (3BB), we obtain (i) discrete three- 
body resonances (3BR), (ii) two-body continuum states 
(2BC) of ^"Li(l+, 2~^)+n, and (iii) three-body continuum 
states (SBC) of ^Li+n+n, which are decomposed from 
the three-body scattering states. Two-body continuum 
states of ^°Li-|-n are expressed by the two straight lines 
whose origins are resonance positions of ^"Li(l+, 2+). 

Since the virtual states of ^°Li(l~, 2~) and of 2n can- 
not be located in CSM due to the limitation of the scaling 
angle 6, the channels of ^''Li(l~, 2~)+n and ^Li-2n com- 
ponents are included in BBC of ^Li-|-n-|-n. 

Using the Green's function, the strength function for 
the operator Ox with rank A is expressed as 



SxiE) = --Im 



3-body scattering (0=0) 
I 



Re(E) 




FIG. 1: Schematic distribution of the energy eigenvalues of 
the ^Li+n+n system with CSM. 



X 6a*.(|') 



(6) 



where Vl'i(^) is the initial wave function of ^^Li. We apply 
the complex scaling to the right hand side of Eq. ©: 



-— Im 

TT 



(7) 



where the complex-scaled Green's function Q^{E, ^, ^') is 
given as 



Q\E,i,i') = {i 

" i E-Et 



E - H{e) 



i'), (8) 
^Ql{E,i,i'iQ) 



In this expansion, E^ and *^(|) (*^(|)) are the 
energy eigenvalues and eigenfunctions (bi-orthogonal 
eigenfunctionsp^ |26^. '2?|') of the complex-scaled Hamil- 
tonian H{9), respectively. Therefore, summation and/or 
integration are taken over v of the solutions of H{9) 
including 3BR, 2BC of ^°U{1+ ,2+)+n and SBC of 
^Li-f n-f (There is no bound state except for the ground 
state). 

Inserting the complex-scaled Green's function in 
Eq. |(7J), we obtain the strength function decomposed into 
each component S\^u{E) for the final state v as 

Sx{E) = ^SxAE). (10) 

" ^t\{dly\^t){^t\6i\^tf 



SxAE) = Im 

TT 



E~E?. 



(II) 



It is noted that the total strength Sx{E) is an observ- 
able and positive definite for any energy and indepen- 
dent of 9. On the other hand, the partial strengths 
Sx,u{E) of resonance and continuum components are not 
necessarily positive definite, and in fact show sometimes 
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negative values. A detailed explanation is discussed in 
Refs.fl^ |2^. Since the number of resonances obtained 
in the calculation depends on 6 and jumps at a certain 
9, Sx^^{E) has a discontinuity at some 9 values p^. 

In the present calculation, we solve the eigenvalue 
problem of the complex-scaled Hamiltonian with the 
hybrid-TV model, and the discretized approximation is 
adopted for continuum states. We use 25 Gaussian basis 
functions for one relative motion and set their maximum 
range about 40 fm. We numerically checked the reliabil- 
ity of the discretized representation of continuum states 
and the stability of the calculated strengths by changing 
the parameters of the basis functions. 

In Fig. 121 we show the eigenvalue distributions of 1/2+ , 
3/2+ and 5/2+ states in a complex energy plane at = 28 
degrees. These are dipole excited states (j = 1) of ^^Li 
from the ground state ( J'^ = 3/2^). We obtain all eigen- 
values along three lines of rotated Riemann cuts corre- 
sponding to two 2BC of i°Li(l+,2+)-hn and one SBC 
of ^Li-|-n-|-n as discussed in Fig. ^ There is no dipole 
resonance which is located between the real energy axis 
and the rotated continuum lines of the excited states. In 
other words, we cannot find any resonances with a sharp 
width at least {r/2Er < tan"! 29; 9 = 28°). This result 
means that the dipole strengths are exhausted by con- 
tinuum states of ^^Li. Since the wave function of s-wave 
valence neutrons is spatially extended, dipole resonances 
including s-wave components in -'^^Li may tend to decay 
easily, and so resonances have large decay widths, whose 
poles are located below the rotated continuum states. 
The effect of such resonances is included in the contin- 
uum spectra in this calculation. 

Using the solutions of the continuum spectra of the 
dipole excited states, we calculate the dipole transition 
from the ground state. In Fig. Ol^a), we show the re- 
sults of the dipole strength functions summing up the 
contributions from the 1/2+ , 3/2+ and 5/2+ states for 
the three types of the ground state wave function. The 
energy is measured from the ^Li-|-n-|-n threshold. It is 
found that the strengths show the low energy enhance- 
ment whose height is sensitive to the (lsi/2)^ probability 




of the ^^Li ground state. This enhancement is interpreted 
as a threshold effect coming from the continuum states 
and reflects the halo structure of ^^Li. We compare our 
results to the experimental data of MSU jj| . The position 
of the enhancement almost agrees with the data, but a 
disagreement of the shape is seen in the strength above 
the energy 1 MeV. We also compare our results with the 
calculation (denoted as DR in Fig. Ota)) by Garrido et 

In Figs. ini[b) and Etc), we derive the cross sections of 
the ^"'^Li breakup by multiplying the transition strength 
and the virtual photon number in the equivalent photon 
method[28l|. where the target is Pb in both cases. We 
can see a good agreement with the data of RIKEnQ for 
the P-2 wave function with the (lsi/2)^ probability being 
around 30%. However, the magnitude of the cross section 
is not observed. It would be desirable to determine its 
magnitude. For the data of GSI 6], the P-1 wave function 
with the (lsi/2)^ probability being around 20% gives a 
good result. However, the experimental error bars are 
still large. Further experimental data with high resolu- 
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FIG. 2: Energy eigenvalues of three dipole excited states 
(1/2+, 3/2+, 5/2+) of ^^Li with 61=28 degree in GSM. Squares 
and triangles indicate 2BG of ^°Li(l+)-|-n and ^°Li(2+)+n, 
respectively. Gircles indicate 3BG of ^Li+n+n. 



FIG. 3: Galculated dipole strengths and cross sections of 
^^Li in comparison to the theory (denoted as DR, Fig. 5 (a) 
of Ref. JlCi|) in (a) and the experimental data; (a) Q, (b)|3 
and (c)Q, where we take into account the convolution. 
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tion and statistics would be required. 

In Fig. 01 we show the result of the separation of the El 
transition strength of ^^Li into two-body and three-body 
continuum components for the three types of the (lsi/2)^ 
probability in the ground state. It is found that the 
two-body continuum component of ^°Li(l+)-|-n shows a 
low energy enhancement in each panel, whose peak posi- 
tion is just above the two-body threshold (0.42 MeV) of 
^''Li(l+)-f n. Another two-body continuum component of 
^°Li(2+)-|-n shows a broader structure because a larger 
decay width of the 2+ state than that of the 1+ state in 
^'^Li broadens the strength. The transition strengths into 
the ^''Li(l+,2+)-|-n 2BC correspond to the following two 
kinds of physical situations: First one is that one of the 
valence neutrons of (Opi/2)^ in the ^^Li ground state is 
excited to a continuum state of s- or d-waves by the El 
external field, and a remaining p-orbital valence neutron 
forms resonances of ^''Li with ^Li. This situation is sim- 
ilar to the case of ^He breakup reactions 14] because the 
^He ground state is almost dominated by (0^3/2)^ of two 
valence neutrons. Second one is that one of the valence 
neutrons of (lsi/2)^ in the ^^Li ground state is excited 
to a p-orbit and forms resonances of ^"Li(l+,2+) with 
^Li, and a remaining s-orbital valence neutron becomes 
a continuum state. This is the characteristic situation of 
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FIG. 4: Decomposition of the dipole transition strength of 
^^Li with the three types of the (lsi/2)^ probability in the 
ground state, as described in the text. 



^^Li. The contribution of this process can be seen as the 
differences between the strengths of 2BC of P-1, P-2, P-3 
in Fig.^ The low energy enhancement in two-body con- 
tinuum components is interpreted as a threshold effect of 
2BC of the ^"Li(l+,2+)+(s-wave neutron) channels. 

From Fig. 01 it is also found that the contribution from 
SBC of ^Li-l-n-fn increases as the (lsi/2)^ probability 
increases. The strength distribution of SBC has a peak at 
the low energy region 0.5 MeV), and slowly decreases 
with energies. At energies higher than about 2 MeV, 
the strength of SBC becomes dominant. The difference 
between the SBC and 2BC strength distributions may 
be understood from the level density of the continuum 
states. The level density of SBC is widely distributed, 
because SBC have two degrees of freedom of the relative 
motions, but 2BC has one degree of freedom. 

With the increasing of the (lsi/2)^ probability in the 
^^Li ground state, the magnitude of the strengths of 2BC 
increases, but the shapes do not change much. On the 
other hand, the strength of SBC shows a sharper peak at 
a low energy, and its magnitude increases markedly with 
the (lsi/2)^ probability in comparison to the case of 2BC. 
This result indicates that the structure of SBC depends 
strongly on the halo structure of the ground state. To see 
it more clearly, we evaluate the integrated strengths up to 
5 MeV for the three types of (lsi/2)^ probabilities, which 
exhaust about 65% of the non-energy weighted sum rule 
value (NEW-SRV) as shown in Table HH We also list 
the ratios of the each component of 2BC and SBC to the 
total integrated strength. It is found that the ratios of 
the contribution of SBC increase proportionally with the 
(lsi/2)^ probability. 

These results mean that the components of 2BC and 
SBC, contribute comparably in the Coulomb breakup 
reaction of ^^Li, and that, in particular, SBC depends 
strongly on the halo structure of the ground state. In 
fact, a large contribution of SBC cannot be seen in the 
case of ^He breakup reactions T^]. In the three-body 
model of ^He, the (lsi/2)^ probability of two valence 
neutrons is 2.4% which is much smaller than that of 
-'^^Li. These results indicate that the mechanisms of the 
breakup reactions of ^He and ""^^Li are different, and this 
fact is caused by the different (lsi/2)^ probabilities in 



TABLE II: Integrated strengths in comparison to NEW-SRV 
and the ratios of the each component of the strength for three 
types of the (lsi/2)^ probability. 





P-1 


P-2 


P-3 


Integral [e^fm^] 


1.401 


1.901 


2.439 


NEW-SRV [e^fm^] 


2.215 


2.856 


3.593 




^Li+n+n 


0.392 


0.522 


0.658 


Ratios 


'°Li(l+)+n 


0.306 


0.262 


0.216 




'°Li(2+)-Fn 


0.302 


0.216 


0.126 
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their ground states. 

There may be broad resonances including virtual states 
which are located under the Riemann cuts of 2BC and 
SBC calculated here. However, their effects are consid- 
ered to produce no remarkable structure, although they 
are seen to give a considerable strength. In the present 
calculation, we cannot conclude how much they could 
contribute to the strength. 

In summary, we investigate the three-body Coulomb 
breakup reaction of ^^Li employing the complex scal- 
ing method (CSM) to describe the three-body unbound 
states of the ^Li+n+n system. We decompose the tran- 
sition strengths into every component of the unbound 
states, such as three-body resonances, and two- and 
three-body continuum states, and examine the effects of 
each component on the strength distribution of ^^Li. 

From the results, we cannot find any dipole resonances 
with a sharp width. The observed low energy enhance- 
ment in the dipole strength comes from the continuum 
states of ^^Li, where the ^"^Li+n two-body continuum 
states and the ^Li+n+n three-body continuum ones give 
comparable contributions. This result also means that 
the breakup mechanism of '^-'^Li is different from that of 
^He where the two-body continuum states correspond- 
ing to the sequential breakup process arc dominant in 
the strength. Furthermore, it is found that the (lsi/2)^- 
component in the ^^Li ground state is responsible for the 
increases of the contribution of the three-body continuum 
states and the low energy enhancement in the strength. 



This calculated Coulomb breakup strength distribution 
of ^^Li shows a good agreement with the experimental 
data of RIKEN. 

As the origin of the low energy enhancement in the 
strength, the threshold effect due to the continuum re- 
sponses is considered, which is seen in both components 
of the two- and three-body contimmm states. Here, the 
strength of three-body continuum states include a con- 
tribution from the virtual states in "'^"Li, which we do 
not consider in this study. It is interesting to investigate 
the effect of virtual states on the breakup reaction in the 
future works. 
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